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We study the behavior of charge carriers in graphene in inhomogeneous perpendicular magnetic 
fields. We consider two types of one-dimensional magnetic profiles, uniform in one direction: a 
sequence of A'' magnetic barriers, and a sequence of alternating magnetic barriers and wells. In 
both cases, we compute the transmission coefficient of the magnetic structure by means of the 
transfer matrix formalism, and the associated conductance. In the first case the structure becomes 
increasingly transparent upon increasing A*' at fixed total magnetic flux. In the second case we find 
strong wave-vector filtering and resonant effects. We also calculate the band structure of a periodic 
magnetic superlattice, and find a wave- vector-dependent gap around zero-energy. 

PACS numbers: 73.21.-b, 73.63.-b, 75.70.Ak 
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I. INTRODUCTION 

The electronic properties of graphene^ in the pres- 
ence of inhomogeneous perpendicular magnetic fields 
have very recently attracted considerable theoretical 
attention^i'^ii^i^i^i^i^'iSiiii. In graphene the charge carri- 
ers close to the Fermi points K and K' form a relativistic 
gas of chiral massless (Dirac-Weyl) quasiparticlcs with, a 
characteristic conical spectrum. This has far-reaching 
consequences. For example, quasiparticles in graphene 
are able to tunnel through high and wide electrostatic po- 
tential barriers, a phenomenon often referred to as Klein 
tunneling, and related to their chiral nature^^. Moreover, 
in a uniform magnetic field graphene exhibits an uncon- 
ventional half-integer quantum Hall effect^'^, which can 
be understood in terms of the existence, among the rela- 
tivistic Landau levels formed by the quasiparticles, of a 
zero-energy one^^. 

From a theoretical perspective, it is then interesting to 
explore how the Dirac-Weyl (DW) nature of the charge 
carriers affects their behavior in non-uniform magnetic 
fields. Such investigation has been started in Ref,-, and 
here we generalize and expand on it by studying several 
more complex geometries. 

Experimentally, inhomogeneous magnetic profiles on 
submicron scales in ordinary 2DEGs in semiconductor 
heterostructures have been produced in several ways, 
and magnetic barriers with heights up to 1 T have been 
obtained. One approach exploits the fringe field pro- 
duced by ferromagnetic stripes fabricated on top of the 
structure^. Another possibility consists in applying a 
uniform magnetic field to a 2DEG with a step^-?. In 
yet another approach, a film of superconducting material 
with the desired pattern is deposited on top of the struc- 
ture, and a uniform magnetic field is applied. In this 
way, magnetic structures with different geometries have 
been experimentally realized, and their mesoscopic trans- 
port properties have been studied, e.g. transport through 
single magnetic barrier^i^ and superlattices^, magnetic 
edge states close to a magnetic step2£, and magnetically 



confined quantum dots or antidots^i. Correspondingly, 
there exists an extensive theoretical literature, pioneered 
by the works of F.M. Peeters and collaborators^-, which 
elucidates the basic mechanisms underlying the behav- 
iors observed in experiments. 

In principle, the same concepts and technologies can be 
used to create similar magnetic structures in graphene, 
once the graphene sheet is covered by an insulating layer, 
which has recently been demonstrated feasible^^*^. Al- 
though at the time of writing there is yet no published 
experimental work demonstrating magnetic barriers in 
graphene, this should be within reach of present-day tech- 
nology, which provides motivation for the present work. 

In a previous paper— we showed that, in contrast 
to electrostatic barriers, a single magnetic barrier in 
graphene totally reflects an incoming electron, provided 
the electron energy does not exceed a threshold value 
related to the total magnetic fiux through the bar- 
rier. Above this threshold, the transmission coefficient 
strongly depends on the incidence angle^iSi. These ob- 
servations were used to argue that charge carriers in 
graphene can be confined by means of magnetic barri- 
ers, which may thus provide efficient tools to control the 
transport properties in future graphene-based nanode- 
vices. 

Here we focus on more complex multiple barrier con- 
figurations and magnetic superlattices. We consider two 
types of one-dimensional profiles. In the first case the 
magnetic field in the barrier regions is always assumed 
to point upwards, while in the second it points alternat- 
ingly upwards and downwards. We shall see that there 
are sharp differences in the transport properties of the 
two cases. 

The outline of the paper is the following. In Sec. [Til we 
introduce the Dirac-Weyl Hamiltonian for graphene, the 
two types of magnetic profiles we consider in the rest of 
the paper, and the transfer matrix formalism for Dirac- 
Weyl particles. In Sec. IIIII and Sec. IIVI we compute and 
discuss the transmission coefficient separately for the two 
cases. In Sec. |V]we consider a periodic magnetic super- 
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lattice, and determine its band structure. Finally, in Sec. 
IVII we summarize our results and draw our conclusions. 



II. HAMILTONIAN AND TRANSFER MATRIX 

Electrons in clean graphene close to the two Fermi 
points K and K' are described by two decoupled copies 
of the Dirac-Weyl (DW) equation. We shall focus here on 
a single valley and neglect the electron spin ."25 Including 
the perpendicular magnetic field via minimal coupling, 
the DW equation reads 



vpcr ■ [ -ihV + -A ) * = i?^. 



(1) 



where cr = [ax,(Jy) are Pauli matrices acting in sublat- 
tice space, and vp = %y. lO^m/s is the Fermi velocity 
in graphene. In the Landau gauge, A — (0, with 
= dx-^^ the y-component of the momentum is a con- 
stant of motion, and the spinor wavefunction can be writ- 
ten as "^{x,y) — '\l){x)e^^y'^ , whereby Eq. ([T]) is reduced 
to a one-dimensional problem: 



-E -idx - i{ky + A{x)) 

-id^ + i{ky + A{x)) -E 



i! = Q. 



(2) 

Eq. ([2|) is written in dimensionless units: with B the typ- 



ical magnitude of the magnetic field, and £b — y/hcfeB 
the associated magnetic length, we express the vector 
potential A{x) in units of Bis, the energy E in units of 
hvp/£B, and x and ky respectively in units of is and Ig^ . 
The values of local magnetic fields in the barrier struc- 
tures produced by ferromagnetic stripes range up to 1 T, 
with typical values of the order of tenth of Tesla. For 
B « 0.1 T, we find £b « 80nm and hvp/iB « 7meV, 
which set the typical length and energy scales. 
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FIG. 1: Magnetic profile Eq. ([Sj: A'^ magnetic barriers of 
width ds separated by nonmagnetic regions of widtfi do. 

We shall consider two types of magnetic field profiles. 
In the first case, illustrated in Fig. [1] the profile consists 
of a sequence of N magnetic barriers of equal height B 



(assumed positive for definiteness) and width ds, sepa- 
rated by nonmagnetic regions of width c?o- The vector 
potential is then chosen as 
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X e [Xn + dB,Xn+l] 

X € [xjv, cx)]. 



(3) 



where n — 0, . . . , N — I and Xn — n(rfo + ds)- The quan- 
tity NdB is the total magnetic flux through the structure 
per unit length in the y-direction.^^ We shall refer to this 
profile as the multiple barriers case and discuss it in Sec. 
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FIG. 2: Magnetic profile Eq. (gjl: iV magnetic barriers of 
width de separated by magnetic wells of width d^B- 

In the second case, illustrated in Fig. [21 each magnetic 
barrier is followed by a region of width d^B of opposite 
magnetic field. The vector potential is accordingly chosen 
as 
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.,N — 1, Xn = n{dB + d- 



(4) 
and 



B- We shall refer to this profile as the 
alternating barrier-well case, and discuss it in Sec. IIVI 
The parameter D has the meaning of net magnetic flux 
through a cell formed by a barrier and a well. For D = 
this profile can be extended to a periodic magnetic su- 
perlattice, a case considered in Sec. IVl 

With our gauge choice, the value of the vector potential 
on the right of the structure is equal to the total magnetic 
flux $ through it 



$ EE A{x > xn) = 



NdB, 
N{dB 



d-i 



case 1, 
case 2, 



(5) 



which is an important control parameter for the transport 
properties. 

In both cases, the solutions to Eq. ^ can be ob- 
tained by first writing the general solution in each region 
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of constant Bz as linear combination (with complex co- 
efficients) of the two independent elementary solutions, 
and then imposing the continuity of the wavefunction at 
the interfaces between regions of different Bz , to fix the 
complex coefficients. This procedure is most conveniently 
performed in the transfer matrix formalism. Here we di- 
rectly use this approach and refer the reader to Refsi^^ 
for a detailed discussion. 
The transfer matrix 



T = 



Til Ti: 
T21 T2; 



(6) 



relates the wave function on the left side of the magnetic 
structure {x < xq = 0) 



ipix) 



1 

-\-iky 



ikl.: 



1 



(7) 



where fc* = y ^ ^y^ ^'^ ^^^'^ wave function on the right 
side {x > xn) 



(8) 




where — — {ky + The coefficients r and 

t are resp. the reflection and transmission amplitudes, 
and we used that with our gauge choices ^ and 
the vector potential vanishes on the left of the magnetic 
structure, and on the right side is equal to <I>. As usual. 



the factor y k^/kl ensures proper normalization of the 
probability current. The relation which expresses the 
continuity of the wave function is then given by^ 



(9) 



Solving Eq. ^ for we get the transmission probability 
T as 




T{E,ky) = \t\ 



H 1 

kl iTnP- 



(10) 



Once T{E, ky) is known, it is straightforward to compute 
the zero-temperature conductance by integrating T over 
one half of the (Fermi) energy surface^i^: 



GiE) = Go 



cos T(i?, £' sin ( 



(11) 



where (j) is the incidence angle (we measure angles with 
respect to the x-direction) , defined by fcj^ = £' sin 0, and 
Go = 2e^ELy/TTh. Ly is the length of the graphene sam- 
ple in the y direction, and Go includes a factor 4 coming 
from the spin and valley degeneracy. 

Before to proceed with the calculations, we can derive 
a simple and general condition for a non- vanishing trans- 
mission. For this purpose it is convenient to parametrize 



the momenta in the leftmost and rightmost regions resp. 
in terms of incidence and emergence angles, (j) and 



ky ~ E sin ( 



kl — Ecostj), 
k{ — E cos 0/ , ky — E sin ( 



(12) 
(13) 



The emergence angle is then fixed by the conservation of 
k-ii'. 



sm( 



sm (Of . 

E 



(14) 



Equation (|14p implies that transmission through the 
structure is only possible if satisfies the condition 



sm (/> + 



< 1. 



(15) 



This condition, already discussed in Ref.^ for the case of 
a single barrier, is in fact completely general and inde- 
pendent of the detailed form of the magnetic field profile. 
It only requires that the magnetic field vanishes outside 
a finite region of space. For |$/£^| > 2, it implies that 
the magnetic structure completely reflects both quasipar- 
ticles and quasiholes. As a consequence of this angular 
threshold, the conductance has an upper bound given by 



Gs{E) = Gn 2- 



e{2\E\-m. 



(16) 



with the Heaviside step function 0. If the vector potential 
profile is monotonous, Gg also coincides with the classi- 
cal conductance, obtained by setting T = 0{1 — \ sin0 -t- 
^/E\). If, however, A{x) is not monotonous, the classical 
conductance is obtained by replacing |$| in Eq. (fTH]) with 
the maximal value of \A\ in the structure, since a classi- 
cal particle is totally reflected as soon as l^lmaa; (rather 
than the total flux |<I>|) exceeds twice the energy. 

Before moving to the next section, we notice, as an 
aside remark, that Eq. ^ can easily be solved in closed 
form for E = Q. The zero-energy spinors are then given 
by 







(17) 
(18) 



These wave functions are admissible if and only if they 
are nornializable, which depends on the sign of ky and 
the behavior of the magnetic field at a; — > ±00. In fact, 
for any A(x), at most one among ^-(_ and ijj^ is admissi- 
ble. If the magnetic field vanishes outside a finite region 
of space, as in our case, one can always choose a gauge 
in which A{x) = on the left of the magnetic region 
and A{x) = $ on its right. It is then straightforward to 
check that for < fcy < — $, the only normalizable solu- 
tion is (fT7|) . whereas for — $ < ky < 0, the normalizable 
solution is psp . In particular, we find that, when the 
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net magnetic flux through the structure vanishes, there 
exist no zero-energy states. This is nicely confirmed by 
the calculation of the spectrum of the periodic magnetic 
superlattice in Sec. |Vl The zero-energy state is a bound 
state localized in the structure. Additional bound states 
of higher energy may also occur-, but we do not further 
investigate this problem here. 



III. MULTIPLE BARRIERS 

In this section we focus on the magnetic profile in Eq. 
In order to compute the transfer matrix T, we need 
the two elementary solutions of the DW equation ^ for 
B, = and the two for B, — 1. 
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FIG. 3: (Color online) Angular dependence of the transmis- 
sion probability at i5 = 1, through N = 1,2,3,10 barriers 
of width ds = 1/A (keeping in this way constant the flux 
$ = Nds ~ 1), and spaced by do ~ 10. 
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FIG. 4: The conductance G/Go at i5 = 1 for several values of 
A'' with d_g = 1/A, such that $ = 1, and do = 10 (black dots 
linked by long-dashed line, which is a guide to the eye only). 
The upper bound (dashed line) corresponds to the classical 
value Gs/Go = 1 (see Eq. (|16p ) while the dotted line is the 
curve given by Eq. (|26|l as a function of N. 



We can then construct the 2x2 matrices Wo{x) and 
Wsix), whose columns are given by the spinor solutions. 
In the nonmagnetic regions we have 

\ E E ^ / 

(19) 

where kj. [x) — ^JE"^ — {ky + A{x))'^. In the regions with 
B^ = I we have 



Wb{x) 



Dpiq) 



p+l 



(20) 



where q ~ \/2{A{x) + ky), p ~ — 1, and Dp{q) is 

the parabolic cylinder function'^^. These matrices play 
the role of partial transfer matrices, and allow us to ex- 
press the condition of continuity of the wave function at 
each interface between the nonmagnetic and the mag- 
netic regions. After straightforward algebra, we get 



(21) 



where 



fn = Wo\^n)WB{Xn)WB^{Xn+dB)Wo{Xn+dB).. (22) 

is the transfer matrix"^" across the {n + 1)"' barrier, and 
we remind that Xn — n{do + cIb)- 

From Eqs. ([22)) . ([2T|) . and pO| we numerically eval- 
uated the transmission probability T for various sets of 
parameters. The results are illustrated in Figs. O [5] and 
El Fig. [3] shows the angular dependence of the transmis- 
sion coefficient at fixed energy for several values of N, 
but keeping constant the magnetic flux <I> through the 
structure. In agreement with the discussion in the pre- 
vious section and Eq. ([TS)) . we observe that the range 
of angles where T 7^ remains the same, 4> < 0, upon 
increasing the number of barriers. At the same time, 
however, the transmission itself is modified, and oscil- 
lations appear, whose number increases with N and for 
larger separations between the barriers. 

More remarkably we find that, rarefying the magnetic 
field by adding more barriers without changing the total 
fiux the transmission probability approaches the clas- 
sical limit, where it is zero or one, depending on whether 
the incidence angle exceeds or not the angular threshold, 
see Fig. O Correspondingly, the conductance as a func- 
tion of N approaches the classical limit, Eq. p^ . see 
Fig. m As expected, the same limit is also approached 
upon increasing the energy, especially for large N. This 
is clearly illustrated in Fig. [5] for the transmission, and 
in Fig. [6] for the conductance: one sees that already for 6 
barriers the classical limit provides a very good approx- 
imation. The classical limit is instead hardly achieved 
changing ds, except when cis is close to the extreme val- 
ues and E/N, see Figs. [I]and[8l 

In conclusion the main result of this Section is that, 
at fixed flux $, the larger is the number of barriers, the 
more transparent is the magnetic structure. This is a 
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FIG. 5: (Color online) Angular dependence of the transmis- 
sion probability, for different values of E, fixing da = 1, 
do = 10, and = 6. 
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FIG. 6: The conductance as a function of the energy for N = 6 
barriers with ds = 1 and do = 10 (solid line). The dashed 
line is the curve given by the classical limit Eq. (|16|l . 



FIG. 7: (Color online) Angular dependence of the transmis- 
sion probability, for different values of ds, fixing E = 2, 
do = 10, and A = 3. 




0.2 0.4 0.6 0.8 1 1.2 

FIG. 8: The conductance as a function of ds for A' = 3 
barriers with E = 2 and do = 10 (solid line). The dashed line 
is the curve given by the classical limit Eq. (|16|l . 



purely quantum mechanical effect, peculiar to magnetic 
barriers. 

Qualitatively this behavior can be explained as follows. 
To calculate the probability for a relativistic particle to 
go through a very thin single barrier, we can simulate the 
profile oi A with a step function with height $/iV in order 
to have the same flux of a magnetic barrier with width 
= ^/N . In this case the transfer matrix is simply 



f = Wo(0")^^>Vo(0+), 



(23) 



where Wo(0~) is given by Eq. with ^ = and 

kx = k^., while in Wo{0'^) we have A — ^/N and con- 
sequently kx — kl. From Eq. (flO)) we get the following 
transmission probability for a single barrier 



^i(</') 



4cos0 COS0/ 61(1 - I sin./) -I- ^/NE\) 
(cos (/) + cos )2 + [$/ {NE)f 



(24) 



where (^/ is defined by Esm(j)f = E'sin^ -I- ^/N . For N 
barriers we can roughly estimate the probability for the 



particle to cross the magnetic structure to be 

TN{<P)^Ti{<j,f B{l-\sin<f> + ^lE\), (25) 

where we have put by hand the global constraint of mo- 
mentum conservation. Using the expression above we 
can then calculate the conductance applying Eq. (|lip . 
To simplify the calculation, in order to have a qualita- 
tive description of the conductance behavior, we further 
approximate cos 4>f — cos (j) in Eq. (j24p , valid for small 
<f>/iV, and replace cosffP' with its average 1/2, getting at 
the end an approximated expression for G which reads 



2N^E^ 



2N^E^ + $2 



N 



(26) 



being $ the total flux for N barriers. In Fig. 0] 
we compare the exact calculation with the approximate 
one given by Eq. (|26p . and find a surprisingly good 
agreement. The small discrepancy has various possi- 
ble sources: the angular dependence of the transmission, 
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which is averaged out in the approximate calculation; the 
finite width of the barriers; the finite separations among 
the barriers, which may let the particles bounce back 
and forth, in this way reducing the transmission. This 
latter effect, therefore, suppresses a bit the conductance 
predicted by Eq. 



IV. 



ALTERNATING MAGNETIC BARRIERS 
AND WELLS 



Next, we consider the magnetic profile in Eq. il- 
lustrated in Fig. [21 In order to construct the transfer 
matrix in this case we need Wb and the partial transfer 
matrix for the regions with = —B, which is given by 



W-b{x) 



Dp+i{~q) 



iV2 



Dp+i{q) 



{p+l)Dp{-q) l^{p+l)Dp{q) 



After some algebra, we then get 



(27) 



■■■fN-2fN-lWZhMm{xN), (28) 

where 

Tn = WZliXn) WB{Xn) Wg^Xn + dg) >V-b(x„ + rffi), 

(29) 

is the transfer matrix2£ across the {n + 1)"^ magnetic 
barrier, and Xn — n{dB + d-B)- Note that Eq. ((29|) 
differs from Eq. (|22p . since now on the right and on the 
left of a magnetic barrier there is a magnetic well rather 
than a nonmagnetic region. 

As in the previous section, the numerical evaluation 
of T is straightforward, and the results for the transmis- 
sion probability and for the conductance are illustrated 
in Figs. El to [111 

Fig. [HI shows the angular dependence of T for a single 
block consisting of a barrier followed by a well of different 
width. The plot emphasizes the very strong wave-vector 
dependence of the transmission, and shows that by tun- 
ing dB and d^B one can achieve very narrow transmit- 
ted beams. This suggests an interesting application of 
this structure as a magnetic filter, where only quasipar- 
ticles incident with an angle within a very small range 
are transmitted. 

More explicitly, in order to select a narrow beam at an 
angle < < -I (at fixed energy), one can choose the 
widths such that 



dB 

d-B 



(1 — sin ( 

2- £2, 



(30) 
(31) 



with £i_2 ^ 1. With this choice, is very close to the 
angular threshold for the first barrier, and the well is 
close to be totally refiecting. The combination of these 
two effects leads to the narrow beam. For — -I < < it 



■" 


\ ' dg=5.8 




dB=2.9 




aB=o.8 


""'"''"^ 


dB=0 



0.5 



0.5 



FIG. 9: (Color online) Angular dependence of the transmis- 
sion probability at i5 = 6 through a structure consisting of 
one magnetic barrier of width da and one magnetic well of 
width d-B- We fix d-B = 11.7 ~ 2E and vary ds from to 
5.8 f» E. 



is enough to flip the magnetic field. These relations hold 
better at small values of |0|. 

In the rest of this section, we focus on structures with 
dB — d-B, i.e. when the total magnetic flux through the 
structure vanishes. 

For several blocks of barriers and wells we observe an 
interesting recursive effect shown in Fig. 1101 There are 
angles at which, upon increasing N , the transmission 
takes at most Nm values, where Nm is the smallest num- 
ber of blocks for which the transmission is perfect (i.e. 
T = 1). Some of these angles are emphasized in Fig. [TO] 
by dashed lines. At those angles, for instance, magnetic 
structures with N equal to integer multiples of 2, 3 and 
5 exhibit perfect transmission, even for small energy. 



N=l ' .. «, fh 
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FIG. 10: (Color online) Angular dependence of the transmis- 
sion probability at iJ = 1 and dB = d-B = I for several 
values of N. The black dashed lines correspond to the angles 
for recursive transmission of multiplicity 2 (at (f) ~ vr/lS), 3 
(at (j) = -7r/6) and 5 (at « -7r/54 and cj> ^ -77r/18). 

These effect can be understood as follows. Suppose 
that, at a given angle, the transmission probability 
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FIG. 11: The conductance G/Go aX E = 1 varying A'^ with FIG. 13: The conductance G/Go as a function of the energy 
dfl = d-B ~ 1. In the inset ds = d-s ~ 0.8 in the same E for A'^ = 6 with ds = d-s = 1 (solid hne). The dashed hne 
range of A*'. is the classical limit, i.e. (2 — ds/E) 6(2E — ds). 
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FIG. 12: (Color online) Angular dependence of the trans- 
mission probability for several values of the energy E, for 
ds = d-B = 1 and A' = 6. 



FIG. 14: (Color online) Angular dependence of the transmis- 
sion probability for several values of ds = d-s with A'^ = 3 
and .B = 1 



through N„i cehs has a resonance and reach the value 
1. Then, for any sequence consisting of a number of cells 
equal to an integer multiple of N^m the transmission is 
1 again. (It is crucial here that, since the magnetic flux 
through each block is zero, the emergence angle always 
coincides with the incidence angle.) At such an angle, 
then, T can only take, upon changing iV, at most Nm val- 
ues. Notice that perfect transmission occurs also for low 
energy of the incident particle, and the angular spreading 
of perfect transmission is also reduced by adding multiple 
blocks. 

This effect in the transmission also reflects in the con- 
ductance. Fig. [Tl] shows that, for a particular set of 
parameters, G oscillates as a function of N with period 
3. However, for a different value of d_B, shown in the 
inset, the period is 2. This unexpectedly strong depen- 
dence of the conductance on adding or removing blocks 
of barriers and wells could be exploited to design a mag- 
netic switch for charge carriers in graphene. Moreover, 
we observe that the angular dependence of the transmis- 



sion is abruptly modified also by changing the energy E 
of the incident particles, see Fig. [T^ or the width of the 
barriers ds. Fig. I14[ where we observe pronounced reso- 
nance effects. As a consequence the conductance exhibits 
a modulated profile as a function of both the energy and 
the barrier's width, as illustrated respectively in Fig. [T3l 
andm 

V. PERIODIC MAGNETIC SUPERLATTICE 

We now focus on the case of a periodic magnetic super- 
lattice. We observe that, if ds = d-B, the profile ^ can 
be extended to a periodic profile, illustrated in Fig. (THl 
where the elementary unit is given by the block formed 
by a barrier and a well. Imposing periodic boundary con- 
ditions on the wavefunction after a length L — xn — N£ 
(where £ ~ 2^^), i.e. iP{xq) — iP{xn), and defining the 
matrix 

n = WB\0)W-BiO)WZl{dB)WB{dB), (32) 
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FIG. 15: The conductance G/Go as a function of ds = d-s 
for A*' = 3 with E — 1. The dashed hne is the classical limit. 
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FIG. 16: Periodic superlattice of alternating magnetic barri- 
ers and wells with ds = d-s- 



standard calculations^^ lead to the quantization condi- 
tion for the energy: 

2 coa{K^e) ^ Tr fl. (33) 

At fixed ky, Eq. (|33p gives the energy as function of 
the Bloch momentum = Notice that Kr^ is re- 

lated to the periodicity of the structure, and parametrizes 
the spectrum. It should not be confused with the x- 
component of the momentum used in the previous 
sections. Fig. [TTl illustrates the first two bands as a func- 
tion of Kx for two values of ky. Fig. [T8]shows the contour 
plot for Tifl as function of E and ky. We find two inter- 
esting main features. First, around zero-energy there is a 
gap, whose width decreases for larger values of j/cyj. This 
is in agreement with the fact that for a magnetic profile 
with zero total fiux there exist no zero-energy states. 

Second, for some values of ky, the group velocity 
Vy^§^ diverges (see Fig. [TH]close to \ky\ > \E\ « 0.3 ). 
The property of superluminal velocity has already been 
observed for massless Klein-Gordon bosons in a periodic 
scalar potential22.. However, to our knowledge, this is 
the first example in which such property is observed for 



massless Dirac-Weyl fermions in a periodic vector poten- 
tial. 



pa 
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FIG. 17: The spectrum for the periodic superlattice in Fig. 
[in] with ds = d-B — 1 Sit ky = (full points) and ky = 1 
(cross poins). is the Bloch momentum. 
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FIG. 18: (Color online) The allowed spectrum, jTrSl| < 2, 
varying E and ky, at ds — d-B ~ 1. The contour lines 
correspond to the values of Tr fl in the interval [—2, 2] at steps 
of 0.5, increasing from blue (inner dark gray) to red (outer 
dark gray). 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have studied the transmission of 
charge carriers in graphene through complex magnetic 
structures consisting of several magnetic barriers and 
wells, and the related transport properties. 



9 



We focussed on two different types of magnetic profiles. 
In the case of a sequence of magnetic barriers, we have 
found that the transparency of the structure is enhanced 
when the same total magnetic flux is distributed over an 
increasing number TV of barriers. The transmission prob- 
ability and the conductance then approach the classical 
limit for large N, see in particular Figs. [3]and|31 

The behavior of alternating barriers and wells turns 
out to be even more interesting. We have shown that a 
single unit consisting of a barrier and a well of suitable 
widths can be used as a very efficient wave- vector filter 
for Dirac-Weyl quasiparticles, see Fig. [S) With several 
blocks we have observed strong resonant effects, such that 
at given angles one gets narrow beams perfectly trans- 
mitted even for low energy of the incident quasiparticles, 
see Fig. (TU] As a result, the conductance is drastically 



modified by adding or removing blocks, see Fig. |4l This 
suggests possible applications as magnetic switches for 
charge carriers in graphene. 

We hope that our paper will further stimulate exper- 
imental work on the rich physics of magnetic structures 
in graphene. 
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